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ABSTRACT 

Context. The dark cloud Lynds 723 (L723) is a low-mass star-forming region where one of the few known cases of a quadrupolar 
CO outflow has been reported. Two recent works have found that the radio continuum source VLA 2, towards the centre of the CO 
outflow, is actually a multiple system of young stellar objects (YSOs). Several line-emission nebulae that lie projected on the east- 
west CO outflow were detected in narrow-band Ha and [S n] images. The spectra of the knots are characteristic of shock-excited gas 
(Herbig-Haro spectra), with supersonic blueshifted velocities, which suggests an optical outflow also powered by the VLA 2 YSO 
system of L723. 

Aims. Our aim is to study L723 in the near-infrared and look for line-emission nebulae associated with the optical and CO outflows. 
Methods. We imaged a field of ~ 5 ' x 5 ' centred on HH 223, which includes the whole region of the quadrupolar CO outflow with 
narrow-band filters centred on the [Fen] 1.644 jim and H2 2.122 jim lines, together with off-line H c and K c filters. The [Fen] and H2 
line-emission structures were identified after extracting the continuum contribution, if any. Their positions were determined from an 
accurate astrometry of the images. 

Results. The Hi line-emission structures appear distributed over a region of 5!5 (~ 0.5 pc for a distance of 300 pc) at both sides of 
the VLA 2 YSO system, with an S-shape morphology, and are projected onto the east-west CO outflow. Most of them were resolved 
in smaller knotty substructures. The [Fen] emission only appears associated with HH 223. An additional nebular emission from the 
continuum in H c and K c appears associated with HH 223-K1, the structure closest to the VLA 2 YSO system, and could be tracing 
the cavity walls. 

Conclusions. We propose that the H9 structures form part of a large-scale near-infrared outflow, which is also associated with the 
VLA 2 YSO system. The current data do not allow us to discern which of the YSOs of VLA 2 is powering this large scale optical/near- 
infrared outflow. 
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1. Introduction 

Lynds 723 (L723) is an isolate dark cloud located at a distance 
of 300 + 150 pc ( |Goldsmith et al.| [1984) 1 that shows evidence 
of low-mass star formation, and is the site where one of the few 
known cases of a quadrupolar CO outflow (two separate pairs of 
red-blue lobes; Lee et al. , 2002, and references therein) has been 



reported. The 3.6 cm radio continuum source VLA 2 ( Anglada 
|et al.| |1996| l, towards the centre of the CO outflow, harbours the 
source that powers the outflow. VLA 2 is embedded in high- 
density gas traced by NH3, which shows evidence of gas heat- 
ing andJme_br^adenmgjGtwF£F^ |1997|l. Two rec ent works 
by |Carrasco-Gonzalez et al.| ( |2008f ar id|Girart et al.|(|2009|l re- 
port that VLA 2 is a multiple system. Carrasco-Gonzalez et al. 
p008] l detect a system of four (VLA 2A, 2B, 2C and 2D) young 
stellar objects (YSOs) and propose that the morphology of the 
CO outflow is actually the result of the superposition of three in- 

Send offprint requests to: R. Lopez 

* Based on observations made with LIRIS at the 4.2 m Williams 
Herschel Telescope operated at the Observatorio del Roque de los 
Muchachos of the Instituto de Astroffsica de Canarias. Fits files of 
Table 1 are available at CDS via anonymous ftp to cdsarc.u-strasbg.fr 
(130.79.128.5) 



dependent pairs of CO lobes. They also propose that one of the 
YSOs (VLA 2A) is powering the largest, east-west pair of CO 
lobes and the system of emission-line optical nebulae, which is 
reminiscent of Herbig-Haro (HH) objects, first reported by Vrba 
|et al.| ( |1986T ). |Girart efaX] p009) detect at 1.35 mm emission 
from dust, resolved into two components, SMA1 and SMA2, 
which have very similar physical properties. SMA2 seems to 
be in a more evolved stage and is harbouring a multiple low- 
mass protostellar system (VLA 2A, 2B and 2C). They also re- 
port emission from the SiO 5-4 line towards the SMA sources, 
this emission shows an elongated morphology that follows the 
northwest-southeast direction pointed by the larger CO outflow 
lobes, which traces a region of interaction between the dense en- 
velope and the outflow. 

In the optical wavelength range, deep Ha and [S 11] narrow- 
band images of the L723 field were obtained by Lopez et al. 



(2006 ). From these images, the knotty structure of the HH 223 
"linear emission feature" (hereafter, we will refer to it as 
HH 223) detected by Vrba et al. ( 1986| l was resolved into sev- 
eral knots, HH 223-A to -F, and the line-emission nature of other 
nebulae in the field was established. Long-slit spectro scopy cov- 
ering the spectral range 5800-8300 A was made by |L6pez et al.| 
(2009) through HH 223. The spectra obtained are characteristic 
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of shock-excited gas (HH spectra), both for the knots and for the 
low-brightness nebula surrounding the knots. The radial veloci- 
ties derived for the knots are supersonic (blueshifted velocities, 
ranging from -60 to -130 km s~'). The velocities derived for 
the low-brightness nebula are compatible with the ambient gas 
velocity. 

In the near-infrared range, Palacios & Eiroa | ( p"999[ l imaged 
a field of ~ 2 ' centred on VLA 2 in the K band and detected H2 
emission from several nebulae located at both sides of VLA 2. 
However, the field imaged only partially covers L723: it does 
not include the complete region encompassed by the CO outflow 
and, in particular, only the western side of HH 223 was mapped. 
Furthermore, since the continuum was not subtracted from their 
narrow-band H2 image, it is impossible to know whether there 
is continuum contribution to the emission in the nebula close to 
VLA 2. Continuum emission around this position is expected if 
the nebula is tracing a cavity opened by the molecular outflow. 

We imaged the L723 field with narrow-band filters centred 
on the [Fen] A 1.644 /mi and H 2 v = 1-0 S(l) (A 2.122 fjm) 
lines together with the broad-band H and the off-line narrow- 
band H c and K c filters. Our aims were to get a more complete 
picture of the HH 223 outflow in the near-infrared wavelength 
range, to search for counterparts of the optical nebulae, and to 
establish the nature of the near-infrared emission coming from 
the nebulae of the L723 field. The main results are presented in 
this paper. 



2. Observations, data reduction and calibration 

Narrow-band images were obtained through filters centred on 
the [Fe n] and H2 v — 1-0 S(l) lines which is adequate to detect 
the emission from shocked ionized and molecular gas respec- 
tively. Additional images were obtained through the broad-band 
H filter and the narrow-band emission-line free filters H c and 
K c , which are useful to subtract the continuum from emission- 
line filters. Observations were made on 2006 July with the 
instrument LIRIS (Long-Slit Intermediate Resolution Infrared 
Spectrograph; Acosta Pulido et al"| 2003} at the 4.2 m Williams 
Herschel Telescope (WHT) of the Observatorio del Roque de los 
Muchachos (ORM, La Palma, Spain). Details about the observa- 
tions are listed in Table [T] LIRIS is equipped with a Rockwell 
Hawaii 1024 x 1024 HgCdTe array detector. The spatial scale is 
0725 pixel -1 , giving an image field of view (FOV) of 4'.21x4'.21. 
The observing strategy was a 5-point dithering pattern. Given the 
elongated morphology of the target we used a E-W offset three 
times larger than that used along the N-S direction. 

The data were processed with the package lirisdr developed 
by the LIRIS team within the IRAF environment] The reduction 
process includes sky subtraction, flat-fielding, correction of geo- 
metrical distortion, and finally a combination of frames using the 
common "shift-and-add" technique. This final step consists of 
dedithering and co-adding frames taken at different dither points 
to obtain a mosaic for each filter. The resulting mosaic covers a 
FOV of ~ 5 arcmin 2 (see Fig. |T 1. Note that narrow-band images 
taken with LIRIS are commonly affected by fringing. An effi- 
cient method to correct this effect is to create a superflat image 
by a median combination of all frames corresponding to a filter 
and rejecting pixels that contain bright stars. 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



Continuum-subtracted images of the emission in the [Fen] 
and H2 lines were obtained by removing the off-line emission 
using the images acquired through the H c and K c filter respec- 
tively. For each pair of H and K images, a flux scaling factor was 
derived by comparing the counts of several stars in both images 
and then subtracting the registered and scaled images. In some 
cases the images with best seeing were degraded by convolution 
with a Gaussian in order to match the width of the point spread 
function (PSF). 

The astrometric calibration of each final image was made 
with the aim of properly comparing the structures at near- 
infrared and optical wavelengths, and with the position of the 
radio continuum sources found in the field. To do that, the near- 
infrared images were registered using the coordinates from the 
2MASS All Sky Catalogue of ten field stars that are well dis- 
tributed over the observed field. The rms of the transformation 
was 0704 in both coordinates. 

In order to estimate the brightness of the detected emission- 
line structures we determined a photometric zero point for each 
filter. For this purpose we obtained instrumental magnitudes in 
each of the final images byusing the image-analysis tools of 
the GAIA interface facilitj[j Firstly, stars in the images were 
detected with the task object detection and then the zero-point 
was determined through a comparison with the magnitudes from 
the 2MASS catalogue. The accuracy reached was better than 0.1 
mag in all the filters. The limiting magnitudes (at a 3<x level) 
found for our images were 20.5 mag arcsec -2 in H c and 18.5 
mag arcsec -2 in K c . Conversion from magnitude to flux density 
was made by using the flux of Vega, which was properly in- 
terpolated taking into account the narrow-band filter response. 
Then, the flux was measured for each substructure that could be 
isolated within the nebular emission features. We defined an el- 
lipse encircling the emission of the substructure to include signal 
down to a ~ 3<x limit by means of the Aperture photometry task 
in GAIA (see Figs . |2] [3} and|4]). The background emission was 
extracted from an ellipse with the same area but encircling a re- 
gion close to the substructure, which was supposedly free from 
line emission. For close substructures, which are difficult to iso- 
late, we preferred to enclose them all together within a single 
ellipse and evaluate the total flux inside. Fluxes are reported in 
Tables|2]and[3] 



3. Results 

3.1. Outflow morphology in H 2 v = 1-0 

Figure [T]displays an image of the L723 field through the H2 v — 
1-0 S(l) narrow-band filter that shows the large-scale emission 
of the HH 223 outflow. As can be seen from the figure, the emis- 
sion spreads over a length of ~ 5'5 (~ 0.5 pc for a distance of 
300 pc) from the southeast to the northwest. It consists of several 
groups of nebulae of different sizes and complexity, distributed 
following a S-shaped pattern that lies projected onto the lobes 
of the east-west CO outflow. The morphology outlined at large 
scale is then suggestive of a parsec-scale H2 outflow. 

The large-scale morphology of the H2 emissions agress well 
with that of the atomic (Ha, [S n]) emissions in the loci where 
both optical and near-infrared emissions are detected. We found 
the H 2 counterparts of the Ha nebulae HH 223-SE1, -SE2, -Nl , 
-N2, and of the largest emission, HH 223 ( |L6pez et al. 2006). 
Additional H2 emission nebulae without visible counterpart ap- 



2 GAIA is a derivative of the Skycat catalogue and image display 
tool, developed as part of the VLT project at ESO. 
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Filter 


AJAA 
W/(A) 


Date 
2006/07/ 


to/r 
(s) 


*exp 

(s) 


seeing 
(arcsec) 


[Fen] 


1.644/280 


15 


100 


2000 


1.05 






18 


100 


600 


0.78 






20 


60 


1800 


0.70 


H 2 u= 1-0 S(l) 


2.122/320 


15 


100 


2000 


1.03 






18 


10 


200 


0.65 






20 


30 


1800 


0.78 


H 


1.625/3100 


15 


5 


500 


1.13 


K, 


2.270/350 


15 


100 


2000 


1.05 






18 


20 


200 


0.75 


H c 


1.570/260 


18 


50 


500 


0.85 



■ 



SB 



Fig. 1. The L723 field imaged with LIRIS through the H2 2.122 fim line filter (continuum is not subtracted). The emission structures 
of the HH 223 outflow have been labeled. North is up and East is to the left. 




Fig. 2. Close-up of the continuum-subtracted H2 image showing 
the structure of the HH 223 emission at the centre of the L723 
field. The emission substructures referred to in Col. 1 of Table|2] 
have been labeled. The white ellipses enclose the regions which 
correspond to the emission-line fluxes referred to in Cols. 3 and 
4 of Table |2] 



pear towards the northwest of HH 223. Some of them were previ- 
ously detected by Palacios & Eiroa ( 1999) (e. g. labeled in their 



work as Kl, K2 and H2-NW). Our new images give a wider cov- 
erage of the HH 223 near-infrared outflow. They include new 
faint emissions (e. g. HH223-NW0, towards the northwest of 



the field) and the filamentary emission K' reported by Hodapp 
( |1994| l, which has a Ha counterpart (HH223-NW2) weaker and 
less extended. In addition, we were able to distinguish the emis- 
sion contributions from the continuum and from the H2 line. The 
continuum-subtracted H2 image shows several emission spots 
that appear to be extended nebulae surrounding embedded stars, 
whose physical association with the stars cannot be established 
from the present data. In contrast, most of the filamentary emis- 
sion features come from "true" line emission that appear com- 
plex and knotty, with several condensations embedded inside a 
more diffuse nebular emission. 

The "true" H2 line emission nebulae along the HH 223 near- 
infrared outflow were identified from the continuum-subtracted 
H2 image, and those with a spatially-extended structure have 
been resolved in smaller substructures. In this work, the H2 neb- 
ulae have been labeled according to the nomenclature of Palacios 
|& Eiroa"| ( |1999| l, but we expanded it with the newly found 
near-infrared features. There is an exception concerning to the 
emission of ~ 30" length towards the centre of the field (i.e. 
HH 223), where we keep the nomenclature established in Lopez 
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Table 2. Positions and photometry of the H2 features detected in the L723 field 



Peak Position Photometry 





^2000 


^2000 


Integrated area 


H2 line flux 




H2 feature 


(h m s) 


(° ' ") 


(arcsec 2 ) 


(10~ 15 erg cm" 2 s~') 


Notes 


SE1 (V83) region 










(1),(4) 


a 


19 18 06.63 


+ 19 11 24.6 








b 


19 18 06.47 


+ 19 11 23.0 








c 


19 18 06.65 


+ 19 11 19.3 








d 


19 18 06.67 


+ 19 11 17.1 








e 


19 18 06.61 


+ 19 11 15.4 








f 


19 18 06.44 


+ 19 11 14.0 








SE2 


19 18 05.48 


+ 19 11 22.1 






(D,(4) 


Nl 


19 17 59.04 


+ 19 11 42.3 


3.77 


0.51+0.08 


(1),(5) 


HH 223 










(2),(6) 


An 


19 17 57.96 


+ 19 11 52.5 






(3) 


A 


19 17 57.96 


+ 19 11 50.1 


7.30 


4.16±0.13 


(2) 


Cn 


19 17 57.55 


+ 19 11 54.8 






(3) 


Ce 


19 17 57.64 


+ 19 11 51.1 


4.68 


6.85+0.10 


(3) 


C 


19 17 57.51 


+ 19 11 52.8 


3.14 


5.58+0.07 


(2) 


Cs 


19 17 57.42 


+ 19 11 50.6 


3.74 


2.78+0.10 


(3) 


E 


19 17 57.04 


+ 19 11 51.5 


4.64 


4.04+0.12 


(2) 


FO 


19 17 56.98 


+ 19 1148.0 


3.94 


0.85+0.07 


(2) 


Fla 


19 17 56.78 


+ 19 11 49.3 






(2) 


Fib 


19 17 56.67 


+ 19 11 50.3 


2.69 


1.26+0.08 


(2) 


H2-N2 






10.32 


0.59+0.15 


(1),(7) 


a 


19 17 55.87 


+ 19 11 52.3 








b 


19 17 55.78 


+ 19 11 52.6 








c 


19 17 55.66 


+ 19 11 53.8 








d 


19 17 55.55 


+ 19 11 54.2 








e 


19 17 55.42 


+ 19 11 54.6 






(3) 


K2 


19 17 54.10 


+ 19 12 13.6 


7.15 


1.54+0.11 


(3) 


Kl 


19 17 53.96 


+ 19 12 18.1 


2.81 


0.53+0.07 




H2-NW 










(3) 


a 


19 17 51.17 


+ 19 12 42.0 


4.24 


12.6+0.10 




b 


19 17 51.30 


+ 19 12 39.6 


4.23 


14.2+0.08 




c 


19 17 51.51 


+ 19 12 36.3 








H2-NW0 


19 17 50.68 


+ 19 12 49.1 






(8) 


H2-NW2 










(9) 


a 


19 1748.11 


+ 19 13 11.4 








b 


19 17 47.56 


+ 19 13 14.4 








c 


19 17 47.09 


+ 19 13 14.8 


7.11 


43.4+0.16 




d 


19 17 46.80 


+ 19 13 13.1 


11.51 


37.3+0.19 




e 


19 17 46.20 


+ 19 13 10.3 


8.97 


33.8+0.21 




f 


19 17 46.11 


+ 19 13 15.0 


11.75 


17.3+0.22 


(10) 


g 


19 17 45.88 


+ 19 13 13.6 









, , Notes: 

(1) With counterpart in Ha; (2) with counterpart in Ha and [S n]; (3) without optical counterpart; (4) position of the apex; (5) position of the peak intensity; (6) this emission has been 
labeled as its optical counterpart; the region between knots C and F corresponds to the H2 emission labeled as H2-SE by Palacios & Eiroa 1 1999 1; (7) knots b to e have been included 
in the flux value; (8) first detection in this work; (9) weak Ha counterpart; detected in K' by Hodapp (1994) ; (10) knots f and g have been included in the flux value. 



|et al.| ( [2006] > for the Ha and [S n] knots. Similar to what we did 
for the Ha and [S n] line-emission features, we added the identi- 
fication label HH223- to all the H2 nebulae found projected onto 
the east-west CO outflow, since they most probably form part of 
the same near-infrared/optical outflow. 

The H2 emission structures identified from our image are 
listed in Table [2] (Col. 1), and their positions are given in Cols. 2 
and 3. The H2 line flux of the nebulae, integrated over the area of 
Col. 4, is given in Col. 5. Figure|2]displays a close-up of HH 223, 
showing the morphology of the H2 emission with detail. Figure 
[3] shows close-up of several extended nebulae (those containing 
smaller substructures) along the HH 223 near-infrared outflow. 

We will briefly describe the main emission nebulae below, 
beginning with HH 223. The rest of H2 emission nebulae are 
listed from southeast to northwest. 

HH 223: The knotty, undulating structure of ~ 22" in length 
found towards the centre of the field, ~ 1' southeast of the ra- 
dio continuum system VLA 2. The H2 feature spatially coin- 



cides with the Ha "linear emission feature" reported by Vrba 
|et al.| ( |T986] >, which corresponds to HH 223 in the Reipurth 
Catalogue of Herbig-Haro objects (Reipurth 1994 1. Part of 
this structure (from HH 223-C to HH 223-F) corresponds to 
the H2-SE emission of |Palacios & Eiroa | ( fl999"l l. The field 



imaged by these authors covers up to the eastern edge of 
HH 223-C, so that HH 223-A was not mapped. Since the 
overall morphology and size of HH 223 is similar for the 
optical and near-infrared emissions and to facilitate the com- 
parison between them, we adopt the nomenclature for the 
knots of |L6pez et al.|(|2006|l, in stead of using the H2-SE la- 
bel given in |Palacios & Eiroa ( 1999| > for the part of the H2 
emission mapped by them. In spite of the similarity between 
the optical and near-infrared emissions, a more detailed com- 
parison shows differences that will be discussed later. 

HH 223-SE1 -SE2: A group of H2 nebulae ~ 2' southeast of 
HH 223-A, the counterparts of the Ha nebulae HH 223- SE1 
(also associated with the V83 nebula of Vrba et al. 1986]) and 
HH 223-SE2, reported by |L6pez et al.| ( |2006| ). HH 223-SE1, 
unresolved in Ha, consists of several arc-shaped structures 
of different sizes (SEl-a to -f) in the H2 image. HH 223-SE2 
appears in H2 as a faint, bow-shaped structure (labeled SE2 
in Tableland Fig. [3}. 

HH 223-N1: Faint, filamentary S-shaped emission that extends 
over ~ 11" in length, ~ 15" southeast of HH 223-A. Only 
the brightest eastern part of it, of ~ 5" in length was barely 
detected in Ha. 
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Peak Position 




Photometry 




Q'2000 


£2000 


Integrated area [Fe 11] line flux 


[Fe n] feature 


(h m s) 


C ' ") 


(arcsec ) 


(10~ 15 ergcirT 2 s -1 ) 


HH 223 (1W2) 










A 


19 17 57.99 


+ 19 11 50.4 


4.68 


5.95 + 0.10 


B 


19 17 57.73 


+ 19 11 49.7 


4.58 


3.39 + 0.10 


C 


19 17 57.45 


+ 19 11 52.9 


5.83 


3.64 ± 0.09 


D 


19 17 57.26 


+ 19 11 52.2 


3.68 


1.16 + 0.09 


E 


19 17 57.05 


+ 19 11 51.8 


3.68 


0.71+0.07 


FO 


19 17 56.96 


+ 19 11 47.2 


2.28 


0.14 + 0.25 


Fla 


19 17 56.77 


+ 19 11 49.6 


1.94 


0.28 ± 0.06 


Fib 


19 17 56.68 


+ 19 11 50.3 


2.18 


0.37 ± 0.07 



Notes: 

(1) The knots were labeled as its optical counterparts. Knot B has not been detected in H 2 ; (2) note the displacements between the [Fen] and 
H2 emission peak positions (see text). 
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Fig. 3. Close-up of the continuum-subtracted H2 image showing 
the other H2 line-emission nebulae of the HH 223 near-infrared 
flow (see Table[2]for identification). The ellipses are the same as 
in Fig. [2] 



HH 223-H2-N2: Faint filamentary H 2 feature ~ 30" northwest 
of HH 223-A and extending over « 15" towards the north- 
west, projected onto the blueshifted CO outflow lobe. This 
emission was barely detected in the image of Pala cios &| 
Eiroa | ( [T999] >. Our H2 image shows several brightness emis- 



sion enhancements (labeled -b to -e in Tableland Fig. [3j. 
No visible counterpart was detected for this emission, except 



for its faint eastern part (H2-N2-a), which coincides with the 
Ha emission HH 223-N2. 

HH 223-K2, -Kl: These two emission nebulae, ~ 6075 (K2) 
and ~ 64" (Kl) northwest of HH 223-A, are the nebulae 
closest to the radio continuum multiple system VLA 2. They 
were first reported by |Palacios & Eiroa | (jT999j as two low- 
brightness H2 patches. K2 appears in our images as a quite 
compact knot of ~ 3" in diameter. Its emission mostly arises 
from the H2 line, as shown in the H2 continuum-subtracted 
frame (see panels in Figs. [2] and [8). The emission of Kl has 
contributions from both continuum and line. In the the H2 
continuum-subtracted frame, Kl shows an elongated shape 
extending ~ 3" along a direction with a position angle 
~ 130° (see panel of Fig. |2]». In contrast its morphology is 
somewhat different in the H2 frame keeping the continuum: 
the Kl emission extends in an emission tail towards its east- 
ern side and, in addition, Kl appears brighter than K2. 

HH 223-H2-NW, -NW0: HH 223-H2-NW is the emission of 
~ 12" length, found ~ V.5 northwest of HH 223-A. Its 
shape follows a smooth undulating pattern, with at least three 
substructures enclosed in a more diffuse emission. North of 
HH 223-H2-NW, an additional filamentary emission of ~ 2" 
length has been detected in our images. It has been labeled 
HH 223-NWO to distinguish it from the rest of the HH 223- 
NW emissions found in the Ha and H2 images, which are 
located few arcmin from HH 223. 

HH 223-H2-NW2: The bright, filamentary H2 emission, with 
several substructures inside, which is seen towards the north- 
western edge of the field. The emission coincides with a 
K' linear nebula reported by Hodapp ( 1994| l. Part of the 
HH 223-H2-NW2 structure is the counterpart of the faint Ha 
nebula HH 223-NW2. 



3.2. [Fen] imaging 

Additional images of the L723 field were obtained in the //-band 
(see Table[T]i. After subtracting the H c from the [Fe n] + contin- 
uum images, significant emission from the [Fe 11] line was only 
detected at the location of HH 223. The structure of the [Fen] 
emission consists of several knotty condensations inside a lower- 
brightness nebular emission (see Fig.Q. In Table [3] we list the 
[Fe 11] condensations of HH 223 (Col. 1), their peak position co- 
ordinates (Cols. 2, 3) and the [Fen] line flux (Col. 5), integrated 
within the area given in Col. 4. As can be seen from the figure, 
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Fig. 4. Close-up of the continuum-subtracted [Fe n] image show- 
ing HH 223. The knotty substructures referred in Table [3] have 
been labeled on the figure. The ellipses are the same as in Figs. 
|2]and|3] 



the [Fe n] emission of HH 223 extends from east to west over 
a similar length than the emissions in the H2 and Ha lines. The 
overall morphologies are also similar. However, some significant 
differences are found when these emissions are compared in de- 
tail. We will discuss this point in the next section. 

Emission from the [Fen] line was not detected out of 
HH 223. Emission from this line has been frequently found 
close to the location of embedded young stellar objects (YSOs), 
and traces the innermost part of the jet that is accelerated near 
the driving source (see, e. g. Reipurth et al.| 2000| i. Hence, one 
would expected to find [Fen] emission associated with Kl, the 
H2 nebula found closest to the radio continuum system VLA 2. 
However, no emission was detected in our images around Kl. 
Likely, the extinction is high enough to prevent the detection, 
provided that the degree of ionization of the outflow is low 
enough to produce only faint [Fe 11] emission. 

It is not reasonable to expect emission from the [Fe n] line 
at regions having Ha emission and no counterpart in the [S n] 
lines (like HH 223-SE or HH 223-NW2). This is because of the 
1 .644 fim [Fe 11] and [S 11] lines arise from similar excited atomic 
levels, and both species have similar ionization potentials. 



4. Discussion 

4. 1 . Comparing the optical and near-infrared emissions of 
the HH 223 outflow 

It is interesting to compare the morphology of the HH 223 out- 
flow emission for different spectral ranges. This can be done 
from the available images of the the L723 field that include the 
outflow. 

At optical wavelengths, we use the deep narrow-band Ha 
and [S n] images acquired with the 2.6-m Nordic Optical 
Telescope (NOT) (see Lope z et aTj |2006| for details of the ob- 
servations). In addition to the near-infrared images in the H and 
K bands discussed in this work, we use an image of the Spitzer 
Infrared Array Camera (IRAC) through its channel 2 (centred at 
4.5 fim). For the typical conditions of low-mass outflows, emis- 
sion from several rotational transitions lines of H2 should be ex- 
pected within the w 4-8 fim wavelength range. It is well-known 
that the IRAC channel 2 is the best suited in detecting shock- 
excited H2 features, in contrast with the IRAC channels 3 (5.8 
fim) and 4 (8.0 fim). This is because channel 2 is the band that is 
less contaminated by extended emission arising from polycyclic 
aromatic hydrocarbon (PAH). The IRAC channel 2 image of the 



L723 field was obtained from the Spitzer Science Archive using 
the Leopard software. The image used here corresponds to the 
post-basic calibrated data (pbcd) obtained in the project POO 139 
From Molecular Cores to Planet (PI N. Evans). 

The original images of the bands we compared have a differ- 
ent spatial scale and FOV. In order to do more meaningful com- 
parisons, all images were transformed to reach the same magni- 
fication. First, we extracted a subimage of the IRAC frame after 
rotating the original pbcd frame by the appropriate angle to be 
oriented in the standard way. Then, we identified six stars that 
are common to all images and are well distributed in the field. 
These stars were used to register the images for converting all of 
them into a common reference system with the iraf tasks geomap 
and geotran. All the final five frames (i. e. Ha, [S n], [Fe 11], H2 
at 2.122 fim and at 4.5 fim) used here to compare the HH 223 
outflow emissions have a spatial scale of 0f.'6 pixel -1 . 

4.1 .1 . Atomic and molecular hydrogen emissions 

Figure [5] displays a composite image of the L723 field, where 
the emission from atomic (Ha) and molecular (H2 2.122 fim 
and 4.5 fim) hydrogen lines appear superposed. The figure re- 
veals that the large-scale morphology of the HH 223 outflow is 
very similar. However, slight differences in several structures can 
be appreciated when the hydrogen emission from each band is 
compared in detail. 

The HH 223-SE1, SE2 and Nl nebulae at the southeastern 
side of the HH 223 outflow are detected in the three bands, and 
they are spatially coincident. However, HH 223-SE2 shows a 
bow-shaped morphology in the two H2 emission bands, while 
this shape is not as clearly outlined in Ha, where the emission 
appears more extended and filamentary. 

To the northwestern side of the outflow, the HH 223-H2- 
NW and H2-NW2 nebulae appear bright and with a very sim- 
ilar morphology in the two H2 emissions, while in Ha H2-NW 
was undetected, and H2-NW2 was barely detected and appears 
shorter. These differences could be due to the extinction, since 
both structures are located within a region that appears highly 
obscured at optical wavelengths. 

The faint filamentary feature HH 223 -N2 detected at 2.122 
fim has neither Ha nor 4.5 fim counterparts. The lack of Ha 
emission a this location should be expected because of the high 
visual extinction. In contrast, we believe that the lack of counter- 
part at 4.5 fim might be better attributed to the weakness of the 
emission, which is below the detection limit of the IRAC frame. 

HH 223, at the centre of the field, shows a morphology and 
a spatial brightness distribution closely coincident for the emis- 
sions in molecular hydrogen at 2.122 fim and 4.5 fim. Some dif- 
ferences appears however when the H2 and Ha emissions are 
compared in detail. One difference is that the low-brightness 
emission surrounding the knots appears more extended in Ha 
than in H2. The knots appear better delimited, with their edges 
sharper in Ha than in H2. Another more relevant difference is the 
lack of counterparts in H2 for knot B, which is very bright in Ha, 
and for knot D, which cannot be resolved as a knotty structure 
within the low-brightness emission, in contrast to what happens 
in Ha. Finally, several faint knotty structures were identified in 
the H2 2. 122 fim line without a counterpart in Ha (i. e. those la- 
beled An, Cs, Ce and F0, this last belonging to the F filament). 
An additional low-brightness filamentary emission in H2 2.122 
fim, without counterpart in Ha, was identified extending ~ 3" 
from knot C to the northeast. The variations of the physical con- 
ditions and, in particular, of the gas kinematics through HH 223 
seem to be mainly responsible for the differences found between 
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Fig. 5. Composite image of the L723 field using a Ha (blue), H2 2.122 fim (green) and IRAC 4.5 fim (red) colour coding to show 
the HH 223 outflow emissions. The positions of the radio continuum sources SMA1 and SMA2 have been marked by plus signs. 
The near-infrared counterpart of the radio source VLA 1 has been marked by a circle. The continua were not subtracted to facilitate 
the alignment of all images using the field stars. 



the atomic and the molecular hydrogen emissions. We will return 
to this point later. 

No Ha counterparts were detected at the positions of the 
HH 223-K1 and K2 nebulae. This should be expected because 
of the high visual extinction in this region, where the radio con- 
tinuum sources are located. In contrast, counterparts at 4.5 fim 
for both HH 223-K1 and K2 are detected in the IRAC frame. 
Note however that at 2.122 fim the emission of HH 223-K2 ap- 
pears rather more extended towards the southeast, away from the 
radio continuum sources, than at 4.5 fim. 

Concerning HH 223-K1, the emission appears significantly 
more extended at 4.5 fim than at 2.122 fim. The emission at 
4.5 /im spreads towards the northwest, beyond its counterpart 
at 2.122 fim, and encompasses the region where the two radio 
continuum sources (SMA1 and SMA2) are located. In contrast, 
as can be seen in Fig. [5] the emission of Kl at 2.122 fim does 
not reach the position of SMA2 (see also Fig. [8]). The extinction 
is probably high enough to prevent its detection in the K band 
around this position. Furthermore, as in the case of the emis- 
sion associated with Kl in the band, its counterpart at 4.5 fim 
could have continuum contribution from heated material of the 
wall cavity opened by the outflow, which can be traced closer to 
the radio continuum sources at this wavelength. 

4.1 .2. [S 11] and [Fe 11] emissions 

The near-infrared [Fen] 1.644 fim emission, when detected, ap- 
pear spatially coincident with the [Sn] A 6716, 31 A emission. 
Both emissions usually have a very similar morphology (see e. g. 
Reipurfh et~aT] |2000| [Davis et al. 



2003). Also, the differences 



between their detailed structures are much smaller than the dif- 
ferences found between the [Fe 11] and H2 emissions. This is be- 
cause the [S 11] and [Fe n] lines arise from atomic levels of similar 
excitation energy and ionization potential. Thus, it is expected 
the emission from both lines to be spatially coincident, tracing 
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Fig. 6. Composite image of HH 223 using a [Sn] 6716,31 A 
(blue) and [Fe 11] 1 .644 fim (red) colour coding (continua were 
not subtracted). North is up and East is to the left. 



the shock-ionized gas component of the outflow. However, the 
relative brightness distribution for both emissions could vary. 
This is not only owing to differential extinction effects, but also 
because of the higher [Fe 11] critical density (^ one order of mag- 
nitude above the [S 11] critical density). Hence, the [Fe 11] emis- 
sion requires a higher pre-shock gas density and/or a faster shock 
velocity to reach an intensity comparable to the [S 11] emission 
intensity. 

Figure [6] shows a close-up of the L723 field including 
HH 223 from knot A to filament F, where the emissions from 
[S n] 67 16,3 1 A (in blue) and [Fe 11] 1 .644 fim (in red) have been 
superposed. The knots detected in the [Fe n] line (see also Fig.Q 
have the corresponding [S 11] knotty counterpart as would be ex- 
pected, although the knots appear more compact in [Fe 11] than in 
[S n]. Moreover, the emission coming from the lower-brightness 
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nebula surrounding the knots is much fainter in [Fe n] than in 
[Sn]. 

Some clues on the spatial distribution and on the relative 
strength of both emissions are obtained from the kinematics and 
physical conditions through HH 223, derived in Lopez et al. 
(2009) from long-slit optical spectroscopy. 

For example, the line-ratio diagnosis shows that knot A (one 
of the brightest in [Fe n]) is the densest, most excited and highly 
ionized knot (indeed, the electron density derived at several po- 
sitions inside knot A fall above the [S n] critical density for col- 
lisional de-excitation). In contrast, the electron densities derived 
in the low-brightness nebula as well as in the fainter HH 223 
knots (D to F), are up to one order of magnitude lower. In ad- 
dition, the excitation and ionization measured in these knots (D 
to F) are visibly lower than in knot A. Hence, the excitation and 
density conditions in knot A could explain its high [Fe n] emis- 
sion as compared with the rest of the knots. This is consistent 
with the results derived by Nisini et al. ( 2002| l and Nisini et 



al. ( 2005) 1 in several HH objects with a more accurate diagnos- 



tic from combined optical/near-infrared spectroscopy. These au- 
thors were able to detect a density stratification from ratios of 
forbidden lines with different critical densities (such as [S n] and 
[Fe n]) as well as variations of the spatial distribution of the ion- 
ization/excitation conditions along the jet emission. In particu- 
lar, they derived that the [Fe n] emission comes from the densest 
and/or highly ionized jet regions. 

On the other hand, note that knot B is bright in [Fen], but 
lacking of apparent emission in H2. Knot C, detected in H2, has 
a [Fen] emission weaker than knot B. Concerning their condi- 
tions, from the optical spectra we derived very similar electron 
densities and excitation/ionization degrees for knots B and C. In 
contrast, knot B is one of the faster HH 223 knots. The radial 
velocity derived for knot B (Vlsr = -125 km s _1 ) is appreciably 
higher (by a factor of two) than for knot C (Vlsr - -55 km s _1 ). 
Thus, we interpret that the discrepancies found between the 
[Fe n] and [S 11] emissions are more likely produced by the kine- 
matics and physical conditions in the knots than by differential 
extinction effects: knots denser and/or faster are stronger emit- 
ters in the [Fen] line than slower and/or more rarefied knots. 
Furthermore, the emission appears spatially more extended in 
[Sn], which gives rise to an extended low -brightness nebula 
around the knots, than in [Fe 11]. This could be explained because 
we are detecting [Fe 11] emission only from the higher density 
material of the knots. Thus, the shape of the overall [Fen] emis- 
sion appears more compact than the [S 11] emission. The lack of 
detection of the H2 line coinciding with the position of the fastest 
knot B also supports this interpretation, since the dissociation of 
the H2 molecule can be more efficient at the knots with higher 
velocity. 

4.1.3. [Fen] and H 2 v = 1-0 S(1) emissions 

Usually, the H2 and [Fe 11] emissions from jets show different 
spatial brightness distributions: the H2 emission is more diffuse 
and extended than the [Fe n] emission (see e. g. Reipurth et al. 



2000). Here we will compare in greater detail the emission from 
both species in the HH 223 outflow. Figure|7]displays a compos- 
ite image, where the emissions from [Fe 11] 1 .644 /im (in green) 
and H2 2.122 fim (in red) appear superposed. 

Several differences between the spatial brightness distribu- 
tion of these two emissions are appreciated. First, the shape of 
some knots is different. This is clearly seen e. g. for knot E. The 
emission of this knot appears more extended in H2 and shows a 
bow-shaped morphology as compared to [Fe 11], where the emis- 



sion is more compact. Second, the relative intensity between the 
knots is also different. Note, e. g. that knots A and C show a sim- 
ilar brightness in H2, while knot C appears appreciably weaker 
in the [Fe 11] line compared to knot A. Even more noticeable is 
the lack of H2 emission at the position of knot B, which is one 
of the brightest in [Fe 11]. Finally, we obtained reliable measure- 
ments of the position offsets between the [Fe 11] and H2 emission 
intensity peaks: In knot A the [Fe 11] emission peaks ~ 077 north- 
east of the H2 emission peak; in knot C the [Fe 11] peak is ~ 1'.' 1 
northwest of the H2 peak (see also Fig. [7] bottom panel). 

Again, we interpret these differences between both emissions 
as mostly caused by the excitation conditions required for each 
line to emit, namely shock velocity and gas density, rather than 
to differential extinction effects. Further support of this interpre- 
tation comes from the anticorrelation that we found between the 
radial velocities of the knots in HH 223, derived from optical 
spectroscopy, and their relative H2 brightness. The radial veloc- 
ity decreases for the sequence of knots B, D, A, E and C, while 
the H2 brightness increases for the sequence of knots C, E, A; 
the emission is very weak in D and is below detection limit in 
knot B where, in addition, the degree of ionization inferred from 
optical line ratios reaches the lowest value measured in all knots. 

As mentioned before, [Fe 11] line emission was only detected 
associated with HH 223. The //-band emission detected around 
the location of Kl, the H2 nebula closest to the radio continuum 
sources, mainly arises from continuum. 

Some of the lower-brightness diffuse H2 emission features 
without [Fe 11] counterpart may trace slower shocks in the ambi- 
ent molecular gas with which the supersonic jet is colliding. This 
was suggested to explain the nature of the H2 emission structures 
lacking of [Fe 11] counterpart that, as in our case, were detected 
along the molecular outflow lobe associated with the jet powered 
by L1551 NE ( |Hayashi & Pyol|2009] ). 

4.2. Near-Infrared environment of the radio continuum 
multiple system VLA 2 

Figure [8] displays a close-up of L723 dark cloud showing the 
field around the low-mass protostellar multiple system (VLA 2A 
to 2D) embedded in the region. The panels of the Fig. [8] display 
the near-infrared emission in the H and K bands obtained from 
the images of this work. 



As mentioned before, Palacios & Eiroa (|1999]> reported 



the detection of two faint H2 nebulae, named Kl and K2, near 
VLA 2. However, our images indicate that the nature of the emis- 
sion is different for Kl and K2. 

Concerning Kl, our images indicate that its emission has 
contributions from both line and continuum. This can be seen 
in the panels of Fig. [3] and its comparison with Fig. [8] The 
morphology of Kl both in the H and K band images includ- 
ing continua consists of a nearly arc-shaped nebula, with two 
(east and west) brightness enhancements of different intensity. 
The emission of the [Fen] 1.644 yum line falls below the de- 
tection limit, as revealed from the continuum-subtracted [Fe 11] 
image. The emission from the H2 2. 122 yum line only contributes 
to the western part of Kl (the side closer to the radio contin- 
uum sources) and appears elongated in the northwest-southeast 
direction (PA =s 140°). Interestingly, |Girart etaLl ( |2009| ) reported 
emission of the SiO 5^1 line towards the position of the radio 
continuum sources. The SiO emission shows an elongated mor- 
phology, reminiscent of a jet, which also follows the northwest- 
southeast direction. Moreover, the H2 emission peaks of Kl and 
K2 are in good coincidence with two SiO emission enhance- 
ments (note that the positions of knots Kl and K2 in Fig. 6 of 
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Fig. 7. Top panel: Composite image, using a H2 2.122 fim (red) 
and [Fen] 1.644 pm (green) colour coding, of a FOV ~ l!5xl!l 
of L723 showing part of the HH 223 outflow (continua were not 
subtracted). Crosses mark the position of the radio continuum 
sources SMA1 and SMA2, close to the H2 structure Kl (see also 
Fig. [8}. Bottom panel: Close-up of the ~ 24'.'3 x 1072 rectangle 
enclosing HH 223. Note the spatial displacements between the 
[Fen] and H2 emissions in the knots. North is up and East is to 
the left. 
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Girart et al. (2009 1 are offset ~ 5" east of the position reported 



in this work). The SiO line is a tracer of shock-excited molecular 
gas. Hence, it seems reasonable to assume that both emissions, 
H2 2.122 pm and the SiO, are tracing the same outflow. The 
continuum emission from the eastern side of Kl may arise from 
heated material of a cavity wall. 

Concerning K2, which lies south of the radio continuum sys- 
tem, the near-infrared emission in the H c and K c (continuous) 
bands falls below the detection limit. In contrast K2 appears 
quite bright, with a magnitude of 17.3, in the H2 narrow-band 
image, and remains without appreciable differences after sub- 
tracting the continuum, both in shape and brightness. Hence, 
K2 emission comes mostly from the 2.122 pm line, which is 
tracing shock-excited gas (note that excitation by fluorescence 
cannot be a likely option because of no ionization sources, like 
early-type stars, are found associated with the L723 dark cloud). 
Interestingly, another clump of SiO emission is detected ~ 6" 
southeast of the K2 intensity peak. This clump is aligned with 
the CO outflow direction. This reinforces the idea that the near- 
infrared (H2 line) and SiO emissions are tracing shock-excited 
gas in a jet, and points towards SMA1 as the likely source pow- 
ering the SiO outflow. 



Fig. 8. Close-up of the L723 field including the K2 and Kl near- 
infrared structures, in the neighbourhood of the radio continuum 
multiple system VLA 2. In the lower-right panel, the crosses 
mark the positions of the sources SMA1 and SMA2 detected at 
1.35 mm by |Girart et al. ( 2009| l. The circles mark the positions 
of the four components of VLA 2 (i. e. VLA 2A to D) detected at 
3.6 cm and 7 mm by Carrasco-Gonzalez et al. (2008). The SiO 
outflow reported by Girart et al. '{ 2009| > has been superposed on 
the image with white contours. The position of the H2 emission 
peak of K2 has been marked with a tilted cross in the other three 
panels. 



4.3. The nature of the radio continuum source VLA 1 



Anglada et al. ( 1991| l detected at 3 .6 cm another radio continuum 
source, VLA 1, located ~ 15" southwest of VLA 2. Later radio 
observations by Anglada lit al.| ( |1996| ) and |Girart et al.| ( | 1997] > led 
to discard VLA 1 as a young stellar object (YSO) associated with 
the CO multipolar outflow. As discussed by Girart et al. ( 2009 ), 
VLA 1 is not associated with the high-density molecular gas, ly- 
ing outside, near the border, of the NH3 structure encompassing 
the radio continuum system VLA 2, and suggest that VLA 1 is a 
line-of-sight radio source unrelated to the outflow. However, the 
nature of the radio source VLA 1 has not been settled up to date. 
Carrasco-Gonzalez et al. ( 2008 1 reanalyse multi epoch, centime- 
tre and millimetre data of L723. Their data are compatible with 
VLA 1 belonging to the L723 cloud and being a radio-emitting, 
optically obscured T-Tauri star. 

Interestingly, we found in all our H- and /f-band images 
a stellar-like, near-infrared counterpart with coordinates a = 
l9 i, 17 ,« 52 s 93; g _ +19 o 12 , y/8, which coincided with i n 0705 

with the position of VLA 1 ( |Carrasco-Gonzalez et al. 2008). 
This counterpart (see Fig. [5} is also detected in the four IRAC 
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channels (3.6, 4.5, 5.8, and 8.0 fim) of the Spitzer images men- 
tioned in Sect. 4. We performed aperture photometry of the 
VLA 1 near-infrared counterpart for all the detected bands and 
found magnitudes of 18.6, 17.0, 15.0, 14.9, 14.8, and 14.8 in 
the He, K c , 3.6, 4.5, 5.8, and 8.0 fim images, respectively (errors 
in magnitude are estimated to be < 0.1 mag). All magnitudes 
were calibrated relative to Vega (see e. g. |Reach et al.| [2005 1 for 
IRAQ. We combined these data to obtain near-infrared colours 
for the VLA 1 counterpart. Following the current classification 
criteria (see e. g. Fang et al. 2009 and references therein), we 



found that the colours of this source are consistent with those of 
an evolved YSO of Class III. In particular, the position of this 
source in the [5.8]-[8.0] vs. [3.6]-[4.5] colour diagram corre- 
sponds to the loci where the Class III YSOs are located, and is 
also compatible within the errors with the loci of the transition 
disc sources. In contrast, the VLA 1 counterpart lies outside the 
diagram region where the background extragalactic sources are 
located. From these data we conclude that VLA 1 is tracing the 
emission of an evolved pre-main sequence star of L723, as pro- 
posed by Carrasco-Gonzalez et al. ( |2008 1. However, VLA 1 is 
likely unrelated to the molecular outflows of L723, because the 
powering sources of these outflows need to be YSOs in an earlier 
evolutionary stage than Class II/III. 

5. Summary and conclusions 

We present near-infrared images of the L723 field, obtained 
with narrow -band filters centred on the [Fen] A 1.644 fim and 
H2 v = 1-0 S(l) lines, together with the off-line filters H c and 
K c . The images cover an area of 4!3 x 4'3, thus include all the 
line-emission nebulae associated with the HH 223 outflow that 
were detected in Ha. The analysis of the near-infrared images 
lead us to the main results summarized below. 

1 . We detected H2 emission in a set of elongated nebular struc- 
tures that appear distributed from the southeast to the north- 
west of the L723 field, extending ~ 5'.5 (~ 0.5 pc for a dis- 
tance of 300 pc), which is reminiscent of a parsec-scale out- 
flow. The H2 structures are found projected on the lobes of 
the larger, east-west CO outflow, with an S-shape morphol- 
ogy. Several H2 structures are new identifications from this 
work. 

2. Additional off-line filter images revealed that there is no 
significant contribution from the continuum to the emission 
of the nebular structures detected in the K band, except in 
HH 223-K1, the structure closest to the position of the radio 
continuum system VLA 2. 

3. Emission from the [Fen] 1.644 fim line was only detected 
in HH 223 at the centre of the imaged field. Additional sig- 
nificant emission in the H band, which mainly arises from 
continuum, was detected in HH 223-K1. 

4. We compared the appearance of the HH 223 outflow in the 
optical (Ha and [S n] lines) and near-infrared ([Fe 11] at 1.644 
fim, and H2 at 2. 122 /mi at 4.5 fim) wavelength ranges, look- 
ing for the nature of the emissions. 

- In general, the HH 223 outflow shows a similar large- 
scale morphology when the emissions from molecular 
(at 2.122 and 4.5 fim) and atomic (Ha) hydrogen lines 
are compared. Some of the discrepancies found among 
the three bands could be attributed to extinction effects 
(e. g. the non detection of Ha counterpart from the near- 
infrared nebulae close to the radio continuum multiple 
system VLA2, and the larger extension of the 4.5 fim 
emission compared with the emission at 2.122 fim, in 



HH 223-K1). Other discrepancies are better explained 
from the physical conditions and kinematics of the re- 
gion where they appear: e. g. the non detection of an H2 
counterpart for the Ha emission at HH 223-B. 

- The [S 11] and [Fe 11] emissions are well coincident for 
all the knots of HH 223 (the only structure of the out- 
flow where near-infrared, line-emission from ionized gas 
was detected). The low-brightness emission surround- 
ing the knots is more extended in [S 11] than in [Fe 11] . 
We attribute this to the lower velocity and density found 
around the knots, which is unable to excite the [Fen] 
emission. 

- Some differences are found when the [Fen] and 
H2 2.122 fim emissions are compared. The differences 
can be better attributed to the different physical condi- 
tions required for these lines to emit (e. g. velocity and 
density) than to be caused by extinction. 

5. We identified the near-infrared counterpart of the radio 
source VLA 1 in the H, K and the four IRAC bands. The po- 
sition of this source in the colour-colour diagrams indicates 
that VLA 1 is tracing the emission of an evolved (Class III) 
YSO. Thus, VLA 1 is likely unrelated to the H 2 outflow. 

6. Finally, we discussed the nature of the near-infrared emission 
for the nebulae closest to the radio continuum sources, which 
is a candidate to power the molecular and optical outflows. 
We propose that the H2 2.122 fim emission in HH 223-K1 
and K2 is tracing shock-excited outflow gas, which also was 
detected in SiO ( |Girart et al.| |2009| ), because the emission 
peaks of the Kl and K2 near-infrared structures are coinci- 
dent with two enhancements of the SiO emission. In addition 
both H2 and SiO emissions appear elongated along a similar 
direction (PA ~ 140°). The extended continuum emission as- 
sociated with HH 223-K1 could be tracing the walls cavity 
opened by the outflow. 

In summary, we propose that the H2 nebular structures de- 
tected in the L723 dark cloud are most probably associated with 
the radio continuum system of protostellar sources found in this 
dark cloud. All these H2 nebular structures could form part of 
a S-shaped outflow that also have an optical counterpart in the 
regions with low visual extinction. 

The knotty structures of the HH 223 outflow could be trac- 
ing internal working surfaces of shocks that originate because 
the gas is ejected at varying speeds and/or with a varying direc- 
tion and faster ejecta overtake slower ones. A variable ejection 
velocity and an undulating jet morphology is expected when the 
exciting outflow source belongs to a binary system. The low- 
brightness, more diffuse line-emission structures could be trac- 
ing slow shocks, excited by the interaction of the CO outflow 
with the accelerated gas of the walls cavity opened by it. The 
emission from the [Fe 11] line could be tracing the densest and 
high-ionized regions in the outflow. 

The current data still do not allow us to discern which of the 
radio continuum component of the VLA 2 system is powering 
the large scale, near-infrared/optical HH 223 outflow. However 
at small scales (i. e. short dynamic time scales), the near-infrared 
emission and the SiO outflow seem to be powered by SMA1. 
The large-scale scenario can be better settled by analysing the 
kinematics and excitation conditions in the line-emission nebu- 
lae and needs to be derived from near-infrared spectroscopy. 
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